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The Holstein Molecular Crystal Model is investigated by a strong coupling perturbative method 
which, unlike the standard Lang-Firsov approach, accounts for retardation effects due to the spread- 
ing of the polaron size. The effective mass is calculated to the second perturbative order in any 
lattice dimensionality for a broad range of (anti)adiabatic regimes and electron-phonon couplings. 
The crossover from a large to a small polaron state is found in all dimensionalities for adiabatic 
and intermediate adiabatic regimes. The phonon dispersion largely smooths such crossover which 
is signalled by polaron mass enhancement and on site localization of the correlation function. The 
notion of self-trapping together with the conditions for the existence of light polarons, mainly in 
two- and three-dimensions, are discussed. By the imaginary time path integral formalism I show 
how non local electron-phonon correlations, due to dispersive phonons, renormalize downwards the 
e-ph coupling justifying the possibility for light and essentially small 2D Holstein polarons. 



1. Introduction 

The interest for phonons and lattice distortions in High 
Temperature Superconductors (HTSc) is today more 
than alive [3, 0]. While the microscopic origin of the 
pairing mechanism in cuprate HTSc has notyet been 
unravelled, evidence has been provided [E 13, IE, @, 0, @| 
that electron-lattice interactions and local lattice fluctua- 
tions are correlated with the onset of the superconducting 
transition [E [IE [H|. Recent investigations [n| would 
suggest a similar role for the lattice also in the newly 
discovered layered pnictide-oxyde quaternary supercon- 
ducting compounds fl3l. [T il [l5f . In fact, the discovery of 
HTSc in copper oxides [16| was motivated by the knowl- 
edge that the Jahn- Teller effect [IE EE EE] is strong in 
Cu 2+ . In nonlinear molecules a lattice distortion lifts the 
degeneracy of the electronic states and lowers the overall 
ground state energy: the energy gain is the Jahn- Teller 
stabilization energy which can be measured optically. As 
vibrational and electronic energies are of the same order, 
the nuclear motion cannot be separated from the elec- 
trons and the combined electron-lattice object becomes 
a mobile Jahn- Teller polaron [2(| which can be described 
in terms of the nonlinear Holstein Hamiltonian 21] . Over 
the last twenty years, the focus on the HTSc has largely 
contributed to trigger the study of the polaron properties 
[221 . 123I Izil . |25| . Earlier and more generally, polarons had 
become a significant branch in condensed matter physics 
[IE HE HE Ell as Landau introduced the concept of an 
electron which can be trapped by digging its own hole in 
an ionic crystal (3(| . A strong coupling of the electron 
to its own lattice deformation implies violation of the 
Migdal theorem [3l[ and polaron collapse of the electron 
band [32[ with the appearance of time retarded interac- 
tions in the system. A great advance in the field came 
after Feynman used the path integrals [33[ to calculate 
energy and effective mass of the Frohlich polaron [34[ by 
a variational method [35[ . In the statistical path integral, 
the quantum mechanical electron motion is represented 
by a real space path x(r) being function of an imaginary 
time t whose range has an upper bound given by the 



inverse temperature. As a key feature of the Feynman 
treatment, the phonon degrees of freedom can be inte- 
grated out exactly (being the action quadratic and linear 
in the oscillator coordinates) but leave a substantial trace 
as a time retarded potential naturally emerges in the ex- 
act integral action. The electron, at any imaginary time, 
interacts with its position at a past time. This self inter- 
action mirrors the fact that the lattice needs some time 
to readjust to the deformation induced by the electron 
motion and follow it. For decades, the Frohlich polaron 
problem has been extensively treated by path integral 
techniques [3ol |37| ] which, remarkably, can be applied for 
any value of the coupling constan t 13811 . A review work 
on the Frohlich polaron is in Ref. [33 ]. 

With regard to small polarons, path integral investiga- 
tions started with the groundbreaking numerical work by 
De Raedt and Lagendijk [IE El[ who derived fundamen- 
tal properties for the Holstein polaron. While a sizable 
electron-phonon coupling is a requisite for polaron for- 
mation j^E HE EiJ] , also the dimensionality and degree of 
adiabaticity of the physical system could essentially de- 
termine the stability and behavior of the polaron states 
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When the characteristic phonon energy is smaller than 
the electronic energy the system is set in the adiabatic 
regime. In materials such as the HTSc, colossal mag- 
netoresistance oxides 60 , [6l| , organic molecular crystals 
[62| . DNA molecules IE Hi] and finite quantum struc- 
tures [IE [sE| , intermediate adiabatic polarons are rele- 
vant as carriers are strongly coupled to high energy op- 
tical phonons. For the HTSc systems, a path integral 
description had been proposed for polaron scattering by 
anharmonic potentials, due to lattice structure instabili- 
ties, as a possible mechanism for the non metallic behav- 
ior of the c-axis resistivity 0, [68| . 

It has been questioned whether small (bi)polarons 
could indeed account for high T c due to their large effec- 
tive mass (69| but, later on, it was recognized that dis- 
persive phonons renormalize the effective coupling in the 
Holstein model yielding much lighter masses [70|. Also 
long range Frohlich e-ph interactions provide a pairing 
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mechanism and bipolaron mass consistent with high T c 
[zElz2,[zS[zH while accounting, among others, for angle- 
resolve d p hotoemission spectroscopy data in cuprates 

htsc izaizl. 

In view of the relevance of thepolaron mass issue, I 
review in this article some work [73] on the Holstein po- 
laron with dispersive phonons which examines the no- 
tion of self-trapping and estimates some polaron prop- 
erties in the parameter space versus lattice dimensional- 
ity. The latter is introduced in the formalism by mod- 
elling the phonon spectrum through a force constant ap- 
proach which weighs the first neighbors intermolecular 
shell. This is done in the spirit of the work by Holstein 
[7^ 1 on the small polaron motion for a one dimensional 
narrow band system. Calculating the polaron hopping 
between localized states due to multiphonon scattering 
processes in perturbation theory, Holstein first pointed 
out that the phonon dispersion was in fact a vital in- 
gredient of the theory whereas a dispersionless spectrum 
would yield a meaningless divergent hopping probabil- 
ity. Time dependent perturbative theory, the hopping 
integral being the perturbation, applies as long as the 
interaction time is shorter than the polaron lifetime on 
a state. Computing the time integral for the hopping 
probability, I had earlier shown [501 ] that such interac- 
tion time is shorter in higher dimensionality thus making 
the time dependent perturbative method more appropri- 
ate in the latter conditions. For a given dimensional- 
ity, the interaction time is also reduced by increasing the 
strength of the intermolecular forces hence, the phonon 
dispersion. Consistently also the crossover temperature 
between band-like motion (at low T) and hopping mo- 
tion (at high T) [7!| shifts upwards, along the T axis, 
in higher dimensional systems with larger coordination 
numbers. Hopping polaron motion is the relevant trans- 
port mechanism in several systems including DNA chains 

It is worth emphasizing that our previous investiga- 
tions on the Holstein polaron and the present paper as- 
sume dimensional effects as driven by the lattice while 
the electron transfer integral is a scalar quantity. In a 
different picture, first proposed by Emin [82| for quasi 
ID solids in the adiabatic limit and recently generalized 
by variational exact diagonalization techniques [831 ] . the 
electron transfer integral is instead a vector whose com- 
ponents are switched on in order to treat anisotropic po- 
laron properties. The system dimensionality is then in- 
creased by tuning the electronic subsystem parameters, 
the 3D case corresponding to an isotropic hopping in- 
tegral, whereas the phonon spectrum is dispersionless. 
While the method to introduce the system dimensional- 
ity in Ref. 83] differs from ours, some trends regarding 
the polaron mass in ID and 3D are apparently at vari- 
ance with ours as it will be pointed out in the following 
discussion. 

Being central to polaron studies, the effective mass be- 
havior ultimately reflects the abovementioned property of 
the polaron problem: when the electron moves through 



the lattice, it induces and drags a lattice deformation 
which however does not follow instantaneously, the re- 
tardation causing a spread in the size of the quasipar- 
ticle. The electron-phonon correlation function provides 
a measure of the deformation around the instantaneous 
position of the electron and may be used to quantify the 
polaron size. Several refined theoretical tools, including 
quantum Monte Carlo calculations (40| . density matrix 
renormalization group [IH , variational methods [H, HE] 
and exact diagonalization techniques [8?], [H[ have been 
applied to the Holstein Hamiltonian to compute such po- 
laron characteristics covering various regimes of electron- 
phonon coupling and adiabatic ratio. While a qualitative 
and often quantitative [84], [H| agreement among several 
numerical methods has emerged especially with regard 
to ground state polaron properties, we feel that analyti- 
cal investigations may provide a useful insight also in the 
specific regime of intermediate adiabaticity and e-ph cou- 
pling for which the leading orders of perturbation theory 
traditionally fail in the weak coupli ng |89l| or become less 
accurate in the strong coupling [90, [9l| . 

Extremely interesting is that range of e-ph couplings 
which are sufficiently strong to allow for polaron forma- 
tion but not too strong to prevent the polaron from mov- 
ing through the lattice. When the relevant energy scales 
for electrons and phonons are not well separated, such 
range of couplings produces a composite particle which 
is not trapped by its lattice deformation. It is this in- 
termediate regime which forms the focus of the present 
work. 

The retardation effect and its consequences are here 
treated by means of a variational analytical method 
based on the Modified Lang Firsov (MLF) transforma- 
tion [92| which considerably improves the standard Lang 
Firsov (LF) [93j ] approach on which strong coupling per- 
turbation theory (SCPT) is based. Section 2 provides the 
generalities of the MLF method applied to the dispersive 
Holstein Hamiltonian. In Section 3, the polaron mass is 
calculated in the (anti) adiabatic regimes covering a broad 
range of cases in parameter space. The spreading of the 
adiabatic polaron size over a few lattice sites is shown 
in Section 4 through computation of the electron-phonon 
correlations. In Section 5, I give further physical motiva- 
tions for the existence of light Holstein polarons: using 
a path-integral method, I show how the electron-phonon 
coupling is renormalized downwards in momentum space 
by nonlocal correlations arising from the dispersive na- 
ture of the phonon spectrum. Some final remarks are in 
Section 6. 



2. Modified Lang-Firsov Method 

The Holstein diatomic molecular model was originally 
cast [2l| in the form of a discrete nonlinear Schrodinger 
equation for electrons whose probability amplitude at a 
molecular site depends on the interatomic vibration coor- 
dinates [HI • The nonlinearities are tuned by the electron- 
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phonon coupling g, whose strength drives the crossover 
between a large and a small polaron for a given value of 
the adiabatic parameter. The polaron radius is measured 
with respect to the lattice constant (95l [96| . When the 
size of the lattice distortion is of the same order of (or 
less than) the lattice constant the polaron has a small 
radius and the discreteness of the lattice must be taken 
into account. 

In second quantization the dimension dependent Hol- 
stein Hamiltonian with dispersive harmonic optical 
phonons reads: 

h= -t cjcj+ff2 + + 12 W 

<ij> i q 

t is the hopping integral and the first sum is over z 
nearest neighbors. c\ and Cj are the real space electron 
creation and annihilation operators at the i-site, rij (= 
c\ci) is the number operator, b\ and bi are the phonons 
creation and annihilation operators. 6 q is the Fourier 

transform of b\ and w q is the frequency of the phonon 
with vector momentum Unlike the Su-Schrieffer- 
Heeger Hamiltonian (97l . |98| a paradigmatic model in 
polymer physics, in Eq. fl} the electron hopping does 
not depend on the relative displacement between adja- 
cent molecular sites hence, the phonon created by b\ is 
locally coupled to the electronic density. 

The LF transformation uses a phonon basis of fixed 
displacements (at the electron residing site) which di- 
agonalizes the Hamiltonian in Eq. (fTJ) in absence of hop- 
ping. T he hopping term is then treated as a perturbation 
70, 90]. However the standard LF approach does not act- 
count for the retardation between the electron and the 
lattice deformations which induces a spread in the size 
of the polaron. Precisely this effect becomes important 
for the intermediate e-ph coupling values which may be 
appropriate for some HTSc. 

The idea underlying the MLF transformation [99( is 
that to consider the displacements of the oscillators at 
different sites around an electron in order to describe 
the retardation effect. For the present case of dispersive 
phonon the MLF transformation, applied to the Hamil- 
tonian in Eq. ([TJ) , reads: 

H = e S He~ S 
S= ^A q n q (&! q -fo q ), 
q 

2 v k 

(2) 

where Ri are the lattice vectors and A q are the varia- 
tional parameters which represent the shifts of the equi- 
librium positions of the oscillators (quantized ion vibra- 
tions) with momentum q. The conventional Lang-Firsov 
transformation is recovered by setting: A q = g/ujq- 



Explicitly, the MLF transformed Holstein Hamiltonian 
in Eq. © is: 

i ij 

xexp[-^^A q 6 q (e^ R »- e ^ R 0] 
q 

x exph-^L £ A q 6 q (e-^ R ' - e"^ R *)] 

+ 12 w q & q & q + 12( g ~ A q w q)"q( 6 -q + & q) 
q q 

(3) 

where the polaron self-energy e p is: 

€ p = TV 12^ g ~ A q^q) A q ( 4 ) 
q 

and the polaronic hopping is: 

* P = *cxp[-l]TA q (l-7 q /z)] 
q 

Tq = 2 12 cosqi ' ( 5 ) 

i=x,y,z 

The coordination number z is twice the system dimen- 
sionality. 

Looking at Eq. ©, it is clear that the unperturbed 
Hamiltonian Hq can be taken as 

H ° = - e P^2 n i + X^qkq^q ' ( 6 ) 

i q 

while the remaining part of the Hamiltonian (H — Hq) 
in the MLF basis is considered as the perturbation part. 
The energy eigenstates of Hq are given by 

\4>i,{ n l}) = c!l >el n qi: n q2>«-q3> :.)ph (7) 

where, i is the electron site and n qi , n q2 , n q3 are the 
phonon occupation numbers in the phonon momentum 
states qi,q2,q3, respectively. The lowest energy eigen- 
state of the unperturbed Hamiltonian has no phonon ex- 
citations, i.e.riq = for all q. The ground state has an 
energy Eq — —e p and is TV-fold degenerate, where N is 
the number of sites in the system. The perturbation lifts 
the degeneracy and to first order in t the ground state 
energy of the 3D- polaron with momentum k is given by 

Eo(k) = -e p -i p 7k (8) 

The second order correction to the ground-state energy 
of the polaron with momentum k is given by 
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x| < {n q },k'|ff-i/ o |k,{0} > | 2 

(9) 

where, in principle, intermediate states having all pos- 
sible phonon numbers contribute to Eq. ([9]). 

By minimizing the zone center ground state energy, the 
variational parameters A q are obtained as: 

Aq = r, TT (10) 

q LUq + Zt p {l - 7q /z) ^ > 

and, by Eq. (|10[) . the one phonon matrix clement be- 
tween the ground state |k = 0, {n q = 0}) and the first 
excited state (l q ,k'| vanishes. Then, the one phonon ex- 
citation process yields no contribution to Eq. ©. By 
Eq. (HH1) one gets the A q 's for the ID, 2D and 3D sys- 
tems and, via Eq. ©, the narrowing of the polaron band 

[loot . 



3. Polaron Mass 

Given the formal background, the polaron mass is cal- 
culated both for the Lang-Firsov and for the Modified 
Lang-Firsov method to the second order in SCPT. Gener- 
ally the second order order correction i ) makes a relevant 
contribution to the ground state energy, ii) does not af- 
fect the bandwidth, Hi) introduces the mass dependence 
on the adiabatic parameter which would be absent in the 
first order. Altogether the polaron landscape introduced 
by the second order SCPT is much more articulated than 
the simple picture suggested by the first order of SCPT 
which nonetheless retains its validity towards the antia- 
diabatic limit. In first order SCPT, ground state energy, 
bandwidth and effective mass appear as equivalent, inter- 
changeable properties to describe the polaron state while 
band narrowing and abrupt mass enhancement are coin- 
cident signatures of small polaron formation in parameter 
space. However, there is no compelling physical reason 
for such coincidence to occur as the bandwidth is in fact 
a zone edge quantity whereas the effective mass is a zone 
center quantity. Thus, by decoupling onset of the band 
narrowing and self trapping of the polaron mass, the sec- 
ond order of SPCT provides the framework for a richer 
polaron structure in momentum space. 

As emphasized in the Introduction, dispersive phonons 
have a relevant role in the Holstein model: a dispersion- 
less spectrum would in fact predict larger polaron band- 
widths in lower dimensionality and yield a divergent site 
jump probability for the small polaron in time dependent 
perturbation theory [78|, Il0l| . Here, a lattice model is 
assumed in which first neighbors molecular sites interact 



through a force constants pair potential. Then, the opti- 
cal phonon spectrum is given in ID, 2D (square lattice) 
and 3D (simple cubic lattice) respectively, by 

2 i \ ct + i 1 r^> — 7. o 

w id (<?) = — + V " + 2a~fcosq + Y 

^L(q) = + ^V« 2 + 2«7.9(q)+7 2 (2 + ^(q)) 

"!d(<0 = + ] J 7 v / « 2 + 2«7j(q)+7 2 (3 + /(q)) 

g(q) = cosq x + cosq y 

h(q) = 2cos(q x - q y ) 

j(q) = cosq x + cosq y + cosq z 

/(q) = 2cos(q x - q y ) + 2cos(q x - q z ) + 2cos(q y - q z )) 

(11) 

where a is the intra-molecular force constant and 7 is 
the inter-molecular first neighbors force constant. M is 
the reduced molecular mass. Thus, the intra- and inter- 
molecular energies are luq = ^J2a/M and u>i — y/j/M 
respectively. Some care should be taken in setting the 
phonon energies as the second order perturbative term 
grows faster than the first order by increasing u>i [70l ]. 
Hence too large dispersions may cause a breakdown of the 
SCPT. In terms of luq, the dimensionless parameter zt/uo 
defines the adiabatic (zt/uo > 1) and the antiadiabatic 
(zt/ojQ < 1) regime. Some other choices are found in 
the literature with t/wo > 1 or tjujQ > 1/4 defining the 
adiabatic regime. As shown below, such discrepancies 
may lead to significantly different interpretations of the 
polaron behavior in parameter space mainly for higher 
dimensions. 

Hereafter I take uiq — lOOmeV and tune t to select a 
set of zt/cuo which sample both antiadiabatic and adia- 
batic regime without reaching the adiabatic limit. The 
dynamics of the lattice is in fact central to our investi- 
gation. A moderate to strong range of e-ph couplings is 
assumed (g/u>o ^ 1) so that the general conditions for 
polaron formation are fulfilled throughout the range of 
adiabatic ratios [22|, [9l| . 

Fig. [T] plots the ratio of the one dimensional polaron 
mass to the bare band mass, against the e-ph cou- 
pling, calculated both in the Lang-Firsov scheme and in 
the Modified Lang-Firsov expression. The intermediate 
regime 2t = ujo is assumed in Fig. [Ua). While at very 
strong couplings the MLF plots converge towards the LF 
predictions, a remarkably different behavior between the 
LF and the MLF mass shows up for moderate g. The LF 
method overestimates the polaron mass for j £ [~ 1 - 2] 
and mostly, it does not capture the rapid mass increase 
found instead in the MLF description. Note that, around 
the crossover, the MLF polaron mass is of order ten times 
the bare band mass choosing uj\ — GOmeV. Large inter- 
molecular energies enhance the phonon spectrum thus 
reducing the effective masses in both methods. In the 
MLF method, large lo\ tend also to smooth the mass be- 
havior in the crossover region. 
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Fig. QJb) presents the case of an adiabatic regime: the 
discrepancies between LF and MLF plots are relevant 
for a broad range of e-ph couplings. The latter is more 
extended for larger intermolecular energies as also seen 
in Fig. [TJa). Mass renormalization is poor in the MLF 
curves up to the crossover which is clearly signalled by 
a sudden although continuous mass enhancement whose 
abruptness is significantly smoothed for the largest values 
of intermolecular energies. 

Fig. [ljc) shows a fully antiadiabatic case in which the 
LF and MLF plots practically overlap throughout the 
whole range of couplings. This confirms that the LF 
method is appropriate in the antiadiabatic limit which is 
essentially free from retardation effects. Unlike the pre- 
vious cases the MLF plots are always smooth versus g 
indicating that no self-trapping event occurs in the an- 
tiadiabatic regime. Here the polaron does not trap as it 
is always small for the whole range of couplings. 

It should be reminded that the concept of self-trapping 
traditionally indicates an abrupt transition between an 
infinite size state at weak e-ph couplings and a finite 
(small) size polaron at strong e-ph couplings. In one di- 
mensio n, accord i ng t o the traditional adiabatic polaron 
theory |l02l . Il03l . Il04l | , the polaron solution is always the 
ground state of the system and no self-trapping occurs. 
Instead, in higher dimensionality a minimum coupling 
strength is required to form finite size polarons, hence 
self-trapped polarons can exist at couplings larger than 
that minimum. 

However, these conclusions have been critically re- 
examined in the recent polaron literature and the same 
notion of polaron size has been questioned due to th e 
complexity of the polaron quasiparticle itself jl05l . Il06l | . 
Certainly, as a shrinking of the polaron size yields a 
weight increase, the polaron mass appears as the most 
reliable indicator of the self-trapping transition. The lat- 
ter, as the results in Fig. [T] suggest, is rather a crossover 
essentially dependent on the degree of adiabaticity of the 
system and crucially shaped by the internal structure of 
the phonon cloud which I have modelled by tuning the in- 
termolecular forces. In this view, self trapping events can 
be found also in the parameter space of ID systems. This 
amounts to say that also finite size polarons can self-trap 
if a sudden, although continuous, change in their effec- 
tive mass occurs for some values of the e-ph couplings in 
some portions of the intermediate/adiabatic regime. The 
continuity of such event follows from the fact that the 
ground state energy is a nalytic function of g for optical 
phonon dispersions [l07j ] hence the possibility of phase 
transitions is ruled out in the Holstein model. Neither at 
finite temperature can phase transitions occu r as the fr ee 
energy is smooth for any phonon dispersion [l08l Il09j | . 

As fluctuations in the lattice displacements around the 
electron site are included in the MLF variational wave- 
function, the calculated polaron mass should not dis- 
play discontinuities by varying the Hamiltonian param- 
eters through the crossover [lldj |. Mathematically the 
crossover points are selected through the simultaneous 
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FIG. 1: Ratio of the one dimensional polaron mass to the 
bare band mass versus e-ph coupling according to the Lang- 
Firsov and the Modified Lang-Firsov methods. The adiabatic 
parameter is set at: (a) the intermediate value, 2t/u>o — 1; 
(b) a fully adiabatic regime, 2t/u)o = 2; (c) an antiadiabatic 
regime, 2t/u>o — 0.25. ujo = WOmeV and ui\ (in units meV) 
are the intramolecular and intermolecular energies of the di- 
atomic molecular chain respectively. 
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occurence of a maximum in the first logarithmic deriva- 
tive and a zero in the second logarithmic derivative of the 
MLF polaron mass with respect to the coupling param- 
eter. Such inflection points, corresponding to the points 
of most rapid increase for m* , are marked in Fig. [2] on 
the plots of the mass ratios computed for a broad range 
of (anti)adiabatic parameters both in one, two and three 
dimensions. 

In ID, see Fig. EJa), the crossover occurs for g values 
between ~ 1.8 — 2.3 and the corresponding self trapped 
polaron masses are of order ^5 — 50 times the bare band 
mass thus suggesting that relatively light small polarons 
can exist in ID molecular solids with high phonon spec- 
trum. The onset of the self-trapping line is set at the 
intermediate value It = ujq and the self-trapped mass 
values grow versus g by increasing the degree of adia- 
baticity. There is no self-trapping in the fully antiadi- 
abatic regime as the electron and the dragged phonon 
cloud form a compact unit also at moderate e-ph cou- 
plings. Then, the mass increase is always smooth in the 
antiadiabatic regime. 

Some significant results are found in 2D as shown in 
Fig. El^b): i) for a given g and adiabaticity ratio, the 
2D mass is lighter than the ID mass and the 2D LF 
limit is attained at a value which is roughly one order of 
magnitude smaller than in ID; ii) the crossover region is 
shifted upwards along the g axis with the self trapping 
events taking place in the range, g ~ 2.2 — 2.6. The 
masses are of order ~ 5—10 times the bare band mass; Hi) 
the line of self-trapping events changes considerably with 
respect to the ID plot. The marked curve is parabolic in 
2D with an extended descending branch starting at the 
intermediate value 4t/ujQ — 1; iv) in the deep adiabatic 
regime, the lattice dimensionality smoothens the mass 
increase versus g. 

This effect is even more evident in 3D, see Fig. [2fc), 
as there are no signs of abrupt mass increase even for 
the largest values of the adiabatic parameter. At the 
crossover, 3D masses are of order ^5 — 10 times the 
bare band mass with the self trapping points lying in 
the range, g ~ 2.5 — 2.9. At very large couplings the 
mass ratio becomes independent of t and converges to- 
wards the LF value. In this region (and for the choice 
uj\ = GOmeV) the 3D Lang-Firsov mass is one order 
of magnitude smaller than the 2D mass. As the coor- 
dination number grows versus dimensionality, large in- 
termolecular forces are more effective in hardening the 
3D phonon spectrum thus leading to lighter 3D polaron 
masses than 2D ones. 

The self trapping transition appears in Fig. [2] as 
smoother in higher dimensions thus contradicting the 
trend found by previous investi gati ons, markedly by ex- 
act diagona lizati on techniques [571 . |83| and variational 
calculations [llll ] . Two reasons may be invoked to resolve 
the discrepancy, the first being more physical and the sec- 
ond more technical. 1) From Eq. (|lip. it can be easily 
seen that: w^[q = (0, 0, 0)] — o^[q = (7r,7r,7r)] = 2dw\. 
Then, for a given uj\, the phonon band is more disper- 
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2: Ratio of the Modified Lang-Firsov polaron 
mass to the bare band mass versus e-ph coupling in 
(a) ID, (b) 2D and (c) 3D. A set of twelve zt/uj 
values ranging from the antiadiabatic to the adiabatic 
regime is considered. From left to right: zt/uoo = 
0.25, 0.5, 0.75, 1.0, 1.25, 1.5, 1.75, 2.0, 2.25, 2.5, 2.75, 3.0. lo = 
lOOmeV. The diamonds mark the occurence of the self- 
trapping event. 
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sive in higher dimensionality d. Hence increasing the 
system dimension corresponds to attribute larger weight 
to the intermolecular forces which ultimately smoothen 
the crossover as mad e evident in Fig. [TJ 2) The works 
in Refs. [57l l83l [TTTI | take t/uio as adiabatic ratio in any 
dimension whereas our calculation spans the same range 
of zt/uo values in any dimension. Thus, for instance, the 
fully adiabatic zt/uio = 3 plots in Fig. [5] would corre- 
spond to an antiadiabatic case in 3D (t/uio = 1/2) and 
a slightly adiabatic case in ID (t/ui p = 3/2) following 
the definition in Refs. [U HU llll| . Consistently, the 
3D antiadiabatic polaron mass is smoother than the ID 
adiabatic polaron mass. Viceversa, assuming the same 
ratio i.e. t/too = 1/2, one should compare the rightmost 
plot in Fig. [2]Jc) with the fourth from left in Fig. EJa): 
the latter is smoother than the former as the ID case is 
now antiadiabatic while the 3D case is fully adiabatic. 

Nonetheless, the findings displayed in Fig.[2]agree with 
density matrix renormalization group [84j . variational 
Hilbert space 15 7J methods and weak coupling pertur- 
bation theory [54j in predicting a lighter mass in higher 
dimensionality. 



4. Electron-Phonon Correlations 

Within the MLF formalism one may also compute 
the electron-phonon correlation functions in the polaron 
ground state. This offers a measure of the polaron size 
as electron and phonons displacements can be taken at 
different neighbors sites. The on site xoi the first neigh- 
bor site Xi an d the second neighbor site \2 correlation 
functions are plotted in Fig. [3] against the e-ph coupling 
both in ID and 2D. For g ~ 1, there is some residual 
quasiparticle weight (not appreciable on this scale) asso- 
ciated with X3 an d Xi which however tend to vanish by 
increasing the coupling. An adiabatic regime is assumed 
to point out the self-trapping event which occurs when 
Xo ~ 1 and Xn ~ for n > 1. This happens for g ~ 2 in 
ID and g ~ 2.5 in 2D for the lowest case of u>i here consid- 
ered. By studying the correlation functions for two values 
of u>i it is seen how the intermolecular forces smooth the 
crossover and extend the polaron size over a larger range 
of g. The case u>i — 60meV, which allows a comparison 
with Fig. [2j is intermediate between the two plots pre- 
sented in Fig. [3j In such case the self trapping appears at 
g ~ 2.3 in ID in fair agreement with the corresponding 
inflection point in the ID effective mass (see fifth curve 
from right in Fig. HI a)). In 2D the polaron size shrinks 
at slightly larger g than in ID and the on site localiza- 
tion occurs at g somewhat larger than the mass inflection 
point. 

Altogether the diamond marked loci displayed in Fig. [2] 
indicate that strong mass renormalization is accompanied 
by on site (or two sites) polaron localization thus iden- 
tifying the self trapping line with the formation of small 
polarons. Before the crossover takes place, for instance 
a t 1 < g < 2 in Fig. [3£a), the ID polaron spreads over a 
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coupling in (a) ID, (b) 2D. An adiabatic regime is considered 
and two values of intermolecular energy are assumed. 



few (two to four) lattice sites. This is the also the range 
of couplings in which the polaron band shows the largest 
deviations [13] from the cosine-like band of standard LF 
perturbation theory [9(| due to the importance of longer 
(than nearest neighbors) hopping processes. The notion 
of intermediate polaron seems to us as the most appro- 
priate for such mobile and relatively light objects. 
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5. Path Integral Method 

The results obtained so far can be put on sound 
physical bases by applying the space-time path integral 
method to the dispersive Holstein Hamiltonian. The 
method permits to incorporate the effect of the electron- 
phonon correlations in a momentum dependent effective 
e-ph coupling. 

The phonons operators in Eq. (fT]) can be generally writ- 
ten in terms of the isotropic displacement field it n as: 

b\ + bt = J? v / 2Afw(q)J]exp( Z q • (R, - R„)K 

q n 

(12) 

Then, the e-ph term in Eq. (TTJ) transforms as follows: 

H e- P h = JLJ2 V2Mw(q) Y, C h U n exp(«q- (Ri - Rn)) 

q i,n 

(13) 

The sum over n spans all n th neighbors of the Ri lat- 
tice site in any dimensionality. It is clear that the phonon 
dispersion introduces e-ph real space correlations which 
would be absent in a dispcrsionless model. Note that 
Eq. (|13p contains the same physics as the e-ph term 
in Eq. ([I]). Fourier transforming the atomic displace- 
ment field and taking the lattice constant |a| = 1, from 
Eq. (TIB)) . I get for a linear chain and a square lattice 
respectively: 

H T Ph = j^Ta E c * c * £ V^MqjexpN' • Ri) 

t q,q' 

x Uq/5d(q / - q) 

n* 

S 1D {q' -q) = l + 2j2 cos(n(q' - q)) 

71=1 

n* 

Saofe' - q) = 1 + 2 Jcos^Aqv) + cos(nAq y ) 

n=l 

n* 

+ 2 j^cos (to Aq x + nAq y ") + cos(mAq x — nAq v ^) 

m,n— 1 

= l'x ~ 1x, Aq y = q' y -q y 

(14) 

While, in principle, the sum over n should cover all the 
N sites in the lattice, the cutoff n* permits to monitor 
the behavior of the coupling term as a function of the 
range of the e-ph correlations. In particular, in ID the 
integer n numbers the neighbors shells up to n* while in 
the square lattice, the n* = 1 term includes the second 
neighbors shell, the sum up to n* = 2 includes the fifth 
neighbors shells, n* — 3 covers the nineth shell and so 
on. Switching off the interatomic forces, w(q) = u>q, one 
would recover from Eq. (| 14[) a local e-ph coupling model 
with Sd = 1- As no approximation has been done at this 
stage Eq. (TT4")) is general. 



A. Semiclassical Holstein Model 

I apply to the Holstein Ha miltonian space-time map- 
ping techniques |112| . Ill3l Ill4 | to write the general path 
integral for an electron particle in a bath of dispersive 
phonons. The m ethod has been used to treat also e-ph 
polymer models [l!5l Ill6j in which the electron hop- 
ping causes a shift in the atomic displacements and, as 
a consequence, the vertex function depends both on the 
electronic and the phononic wave vector. Such shift is 
howe ver absent in the Holstein model as Eq. (fl]) makes 
clear [H3,[ii3. 

I introduce x(r) and y(r') as the electron coordinates 
at the i and j lattice sites respectively, and the electronic 
Hamiltonian (first term in Eq. |l|) transforms into 

iJ e (r,r') = -t( C t(x(r))c(y(r')) + C t(y(r'))c(x(r))) 

(15) 

r and t' are continuous variables (e [0, (3\) in the Mat- 
subara Green's functions formalism with /3 being the in- 
verse temperature hence the electron hops are not con- 
strained to first neighbors sites. After setting r' = 0, 
y(0) = 0, the electron operators are thermally averaged 
over the ground state of the Hamiltonian. As a result, the 
average energy per lattice site, h e (r) = < H c '(t) > /N, 
associated to electron hopping reads (in d dimensions): 

h e (r) = -<(G[-x(r),-r] +G[x(t),t]) 

G[xW , Tl ^/^ rfk . x(T)1? ^=^ 

(16) 

v n are the fermionic Matsubara frequencies and ek is 
the electron dispersion relation. 

Consider now the e-ph term. The spatial e-ph corre- 
lations contained in Eq. (|14[) are mapped onto the time 
axis introducing the r dependence in the displacement 
field: M q — * Uq(r). Assuming periodic atomic particle 
paths: u q (r + /3) = it q (r), the displacement is expanded 
in Np Fourier components: 

N F 

u q (r) = u + 2 f fflu w )q cqs(oj»t) - (gm n ) q sm(uj n T)] 

n=l 

(17) 

with uj n — 2nir/f3. I take a semiclassical version of the 
Holstein Hamiltonian [l!9| . assuming that the phonon 
coordinates in Eq. (fl"4"]) as classical variables interacting 
with quantum mechanical fermion operators. Such ap- 
proximation may affect the thermodynamics of the sys- 
tem as the quantum lattice fluctuations in fact play a 
role mai nly for intermediate values of the e-ph coupling 

Sing. 

Averaging Eq. (fT4")l on the electronic ground state, the 
e-ph energy per lattice site is defined as 
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< HT vh > 

N 



E<^ e 



e—ph 



>c 



< H e d - ph >q= ^72 V^MqJE^'^'^^' - 



q' 



(18) 



Then, on the base of Eq. (ITTD a nd Eq. [fig]). I identify 
the perturbing source current 121] for the Holstein model 
with the t dependent averaged e-ph Hamiltonian term: 



j( T ) = Eiq( r ) 



H {r)^<Hr ph > c 



(19) 



Note that the Holstein source current does not depend 
on the electron path coordinates. The time dependence 
is incorporated only in the atomic displacements. This 
property will allow us to disentangle phonon and electron 
degrees of freedom in the path integral and in the total 
partition function. 

After these premises, one can proceed to write the gen- 
eral path integral for an Holstein electron in a bath of dis- 
persive phonons. Assuming a mixed representation, the 
electron paths are taken in real space while the phonon 
paths are in momentum space. Thus, the electron path 
integral reads: 



< x(/3)|x(0) >= J] < x (/ 3 )l x (°) >c 



< x(/3)|x(0) > q = / Duq(r)exp 



W(q)<(r)) 



x J Dx{ T )exp - jy T ^\r)+h e {r)-U{r)\ 

(20) 

where m is the electron mass. The perturbing cur- 
rent in Eq. (|20|) is integrated over r using Eq. (fl7| and 
Eq. IHl). The result is: 



drjqir) = (3u g d (q) 



gM) = ^V2Ma;(q)5> q ,S d (q' - q) 
q' 



(21) 



where <?d(q) is thus a time averaged e-ph potential. 

The total partition function can be derived from 
Eq. (|20|) by imposing the closure condition both on the 
phonons (Eq. (fTT)) ') and on the electron paths, x(/3) = 
x(0). Using Eq. J2l]), I obtain: 



M 



q J 

+ ulu\{r)) 
x <z> Dxexp 



(3u g d (q) - / dr— (u' q (r) 



(22) 



§ Duq and § Dx are the measures of integration which 
normalize the kinetic ter ms in the phonon field and elec- 
tron actions respectively [l2l| . 

The phonon degrees of freedom in Eq. (j22|) can be in- 
tegrated out analytically yielding: 



Z T = H>(q) x j D *i T ) 



x exp 



P(q 



1 



■ exp 



(sd(q)A 



M 



ff (2n7r) 2 

1 = 1 (2nvr)2 + (/3 w (q))S 



2vrw(q) 5 



(23) 



being X M = y/irh 2 (3/M. 

Eq. (|23|) is the final analytical result from whic h th e 
thermodynamics of the model can be computed [ll8j |. 
The exponential function in P(q) embodies the effect of 
the non local correlations due to the dispersive nature 
of the phonon spectrum. Phonon and electron contri- 
butions to the partition function are decoupled although 
the effective potential <?d(q) carries a dependence on the 
electron density profile in momentum space through the 
function pq. 



B. Electron-Phonon Coupling 

The time (temperature) averaged e-ph potential in 
Eq. (|2Tj) is computed in the case of a linear chain and 
of a square lattice. As an example, I take the elec- 
tron density profile as p q — p Q cos(g) in ID and p q = 
p cos{q x ) cos{q y ) for the 2D system. Since the momen- 
tum integration runs over ^ £ [0,7r/2], p represents in 
both cases the total electron density. More structured 
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FIG. 4: (a) Time averaged e-ph coupling (in units meVA~ 
versus wave vector for a linear chain, n* represents the cutoff 
on the e-ph correlations, po is the electron density. u>o = 
20meV and cut — WmeV. g = 3. The dispersionless e-ph 
coupling is obtained for u/i = 0. (b) Time averaged e-ph 
coupling in two dimensions versus the y component of the 
momentum. The cases n* = 1, n* = 2 and n" — 3 imply that 
the correlation range includes the second, the fifth and the 
nineth neighbors shell, respectively. The input parameters 
are as in (a). 



density functions containing longer range oscillations are 
not expected to change the trend of the results hereafter 
presented. Low energy phonon spectra parameters are 
here assumed, loq = 20meV and u)\ = lOmel^. Setting 
g = 3, I take a strong Holstein coupling although the 
general trend of the results holds for any value of g. 

In Fig. [3Ja), gio{q)/po is plotted for three choices of 
the cutoff n* to emphasize how the potential depends on 
the e-ph interaction range. The constant value of the 
potential obtained for u>i = is also reported on. For 
short correlations (n* = 4) there is a range of wave vec- 
tors in which the effective coupling becomes even larger 
than in the dispersionless case: this is due to the fact 
that the approximation related to a short cutoff is not 
sufficiently accurate. In fact, extending the range of the 
e-ph correlations, the effective potential is progressively 
and substantially reduced for all momenta with respect 
to the dispersionless Holstein model. Numerical conver- 
gence for <?i_d(q) is found at the value n* = 24 which 
corresponds to 48 lattice sites along the chain. Then, the 
ID coupling renormalization is q-dependent and gener- 
ally larger for q ~ tt/2 where the phonon dispersion gen- 
erates stronger e-ph correlations. Exact diagonalization 
techniques [87| had found that multiphonon states with 
longer range hopping processes substantially soften the 
polaron band narrowing with respect to the predictions 
of standard SCPT with the largest deviations occuring 
at the polaron momentum ~ tt/2. In fact strong, non- 
local e-ph correlations should favour hopping on further 
than neighbors sites occupied by phonons belonging to 
the lattice deformation. 

The projections of the two dimensional e-ph potential 
along the y component of the wave vector is shown in 
Fig. [3fb) for three values of the cutoff. For n* = 1 the 
correlation range is extended to the second neighbor shell 
thus including 8 lattice sites. For n* = 2 and n* — 3 the 
normalization is over 24 and 48 lattice sites respectively. 
Then, the three values of n* in 2D span as many lattice 
sites as the three values of n* in ID respectively. I have 
made this choice to normalize consistently the potential 
for the linear chain and the square lattice. There is a 
strong renormalization for the 2D effective potential with 
respect to the dispersionless case for any value of the wave 
vector. Remarkably, the reduction is already evident for 
n* = 1 . By extending the correlation range this tendency 
becomes more pronounced for q y close to the center and 
to the edge of the reduced Brillouin zone. An analogous 
behavior is found by projecting the 2D potential along 
the q x axis. The 2D potential stabilizes by including the 
9th neighbor shell (n* = 3) in the correlation range. 

Then, an increased range for the e-ph correlations leads 
to a very strong reduction of the effective coupling, an ef- 
fect which is qualitatively analogous to that one would 
get by hardening the phonon spectrum. This is the phys- 
ical reason underlying the fact that the Holstein polaron 
mass can be light. 
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6. Conclusions 

According to the traditional notion of small polaron, 
the strong electron coupling to the lattice deformation 
implies a polaron collapse of the electron bandwidth (zt) 
and remarkable mass renormalization. Thus, bipolaronic 
theories have been thought of being inconsistent with su- 
perconductivity at high T c as the latter is inversely pro- 
portional to the bipolaron effective mass. Motivated by 
these issues, I have investigated the conditions for the 
existence of light polarons in the Holstcin model and ex- 
amined the concept of self-trapping versus dimensional- 
ity for a broad range of (anti)adiabatic regimes. The self 
trapping events correspond to the points of most rapide 
increase for the polaron effective mass versus the e-ph 
coupling. A modified version of the Lang-Firsov trans- 
formation accounts for the spreading of the polaron size 
due to retardation effects which are relevant mainly for 
moderate e-ph couplings. Assuming large optical phonon 
frequencies (u>o), it is found that polaron masses display 
an abrupt although continuous increase driven by the e- 
ph coupling in any dimension, included the ID case. Such 
self trapping events occur in the adiabatic and intermedi- 
ate (zt/too ~ 1) regimes provided that the e-ph coupling 
is very strong (g/ujQ > 2). The second order of strong 
coupling perturbation theory permits to decouple two 
fundamental properties of the polaron landscape, namely 
band narrowing and effective mass enhancement. I have 
considerd a range of e-ph couplings (g/wo > 1) in which 
the conditions for polaron formation are fulfilled. For a 
given value of the adiabatic ratio, the onset of the band 
narrowing anticipates, along the g axis, the crossover to 
the heavy polaron state which pins the quasiparticle es- 
sentially on one lattice site. In such buffer zone, between 
band narrowing onset and self trapping, the charge carri- 



ers can be appropriately defined as intermediate polarons. 

In all dimensions, there is room for non trapped inter- 
mediate polarons in the region of the moderate couplings. 
Such polarons spreads over a few lattice sites and their 
real space extension grows with dimensionality. Thus, 
analytical methods can suitably describe some polaron 
properties also in that interesting intermediate window 
of e-ph couplings (i.e. 1 < g/ojQ < 2 for the ID system) 
in which adiabatic polarons are not large but not yet on 
site localized. In the square lattice, at the onset of the 
self trapping state, the polaron over bare electron mass 
ratios are ~ 5 — 10 and even smaller for lower couplings. 

The physical origins for the possibility of light po- 
larons have been analysed treating the dispersive Hol- 
stein Hamiltonian model by the imaginary time path in- 
tegral method. The perturbing source current is the time 
averaged e-ph Hamiltonian term. The partition function 
of the model has been derived integrating out the phonon 
degrees of freedom in semiclassical approximation. Fo- 
cusing on the role of the intermolecular forces, I have 
quantitatively determined the effects of the dispersive 
spectrum on the time-averaged e-ph coupling which in- 
corporates the nonlocality effects. It is found that, in the 
square lattice, the e-ph coupling is strongly renormalized 
downwards and this phenomenon is already pronounced 
by the (minimal) inclusion of the second neighbors in- 
termolecular shell. This may explain why 2D Holstein 
polarons can be light while the real space range of the 
e-ph correlations is relatively short. 
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